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Abstract 

The paper is devoted to the connection between integrabihty of a finite quantum system and 
degeneracies of its energy levels. In particular, we analyze in detail the energy spectra of finite 
Hubbard chains. Heilmann and Lieb demonstrated that in these systems there are crossings 
of levels of the same parameter independent symmetry. We show that this apparent violation 
of the Wigner-von Neumann noncrossing rule follows directly from the existence of nontrivial 
conservation laws and is a characteristic signature of quantum integrabihty. The energy spectra 
of Hubbard chains display many instances of permanent (at all values of the coupling) twofold 
degeneracies that cannot be explained by parameter independent symmetries. We relate these 
degeneracies to the different transformation properties of the conserved currents under spatial 
reflections and the particle-hole transformation and estimate the fraction of doubly degenerate 
states. We also discuss multiply degenerate eigenstates of the Hubbard Hamiltonian. The wave 
functions of many of these states do not depend on the coupling, which suggests the existence 
of an additional parameter independent symmetry. 

1 Introduction 

The close connection between symmetry and degeneracy has been explored since the foundation of 
Quantum Mechanics. Famous examples include degeneracies of spectra in angular momentum in the 
Hydrogen atom (well known as the accidental degeneracy) and the 3d harmonic oscillator [0]. In 
the present paper we focus on the implications of the rich symmetry structure for the spectrum of the 
Id Hubbard model as well as for a general class of quantum integrable systems. We analyze in detail 
how degeneracies and other spectral properties of the Hubbard Hamiltonian reflect the symmetries 
^ , of the model. 

Here we restrict ourselves to the case of Hamiltonians that depend on a single real parameter, 
referred to as a coupling. In this case one can distinguish between two types of degeneracy. The 
first type, often called permanent degeneracy, refers to energy levels which remain degenerate for 
all values of the parameter. The appearance of permanent degeneracies suggests the existence of 
non-commuting symmetry operators (see e.g. [^). 

The second type of degeneracy is a crossing of energy levels that occurs at a particular value of 
the coupling. Multiple level crossings at a certain point indicate a higher symmetry of the system at 
this particular point as compared to other points. For example, the limits of zero or infinite coupling 
frequently have enhanced symmetry. 
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Simple pairwise crossings are also interesting from a different perspective. As is well known, in the 
absence of symmetries energy levels repel, or, equivalently, crossings of levels are prohibited unless 
the states have different symmetry. This noncrossing rule was suggested by Hund in 1927 0, and 
justified by von Neumann and Wigner |Q two years latter. Since then the level repulsion phenomenon 
has been revisited and elaborated upon in various contexts, see for example ^, |^, and references 
therein. Usually textbooks, e.g., ||, ||, |], present a simplified version of the justification due to Teller 

The mathematical validity of the noncrossing rule depends crucially on the interpretation of the 
word symmetry, normally understood as a space (or internal space) symmetry. This was highlighted 
by an interesting example by Heilmann and Lieb in 1971 |Tl|. These authors pointed out that the 
noncrossing rule is apparently violated in the case of the Id Hubbard Hamiltonian for the benzene 
molecule if only those symmetries that do not depend on the coupling constant are taken into account. 
Indeed, Fig. ^ illustrates that a substantial number of crossings can be found even within the subsets 
of levels characterized by the same set of quantum numbers. Based on this fact, Heilmann and Lieb 
concluded that "there must be a natural parameter dependent group" to account for these violations. 

The Id Hubbard model is solvable by Bethe's Ansatz, and usually the applicability of this Ansatz is 
understood to rest upon Quantum Integrability of the model, i.e the existence of an infinite number of 
mutually commuting operators. Such commuting operators are often alternately termed as "conserved 
currents " , or "dynamical conservation laws" , or simply "dynamical symmetries" , and are invariably 
associated with all known quantum integrable models. However, it was not until 1986 when the 



parameter dependent integrals of motion (conserved currents), were identified by Shastry |16, 17, 18 
who constructed a transfer matrix that commuted with the Hubbard Hamiltonian. 
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Figure 1: Energies for (P, cr, S, L) = (3, 1, 1, 0, 1) in units oi U — AT as functions oi u = U /{U — 
AT). Levels marked with crosses are twofold degenerate. They also correspond to {P,a,I^°\ S,L) = 
(3, —1, 1, 0, 1). Upside down the figure shows levels for the same set of quantum numbers with S ^ L. 
***** shows the seven fold degenerate level. 

The Hubbard Hamiltonian to be considered consists of the kinetic term, which allows electrons 
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to hop between the nearest neighbors on a regular polygon (i.e. periodic boundary conditions are 
assumed) and the Coulomb interaction between electrons of opposite spin on the same site. 

H = TJ2J2 (cj,c,+i, + cj+i,c,,) + f/E ^.T - ^ %i - ^ (1-1) 
j=is=u i=i ^ ^ 

where c]^ and Cjs are the creation and annihilation operators for an electron of spin projection s on 
site j and hjg = c^jgCjs is the number operator. The goal of this paper is to study how the interplay 
between parameter dependent integrals of motion and ordinary (parameter independent) symmetries 
is manifested in permanent degeneracies and level crossings in the energy spectrum. The Hamiltonian 
( p. . 1| ) is chosen as a tutorial example of a many body system that has both parameter dependent and 
parameter independent symmetries. 

The paper is organized as follows. Section 2 is devoted to parameter independent symmetries 
of the Hubbard model that are later used in Section 4 to diagonalize the Hamiltonian and assign 
symmetry quantum numbers such as the total momentum, spin, particle-hole symmetry etc. to all 
states. Section 3 discusses the structure and properties of parameter dependent conservation laws. 
Permanent degeneracies are treated in detail in Section 5. In the same section we consider the multiple 
permanent degeneracies and mention spectral properties of the dynamical conservation laws. 

In Section 6 we formulate the noncrossing rule in the framework of adiabatic equations of motion 
for matrix elements of the Hamiltonian and higher currents. This language is useful for understanding 
the suppression of level repulsion in integrable models. Next, based on numerical analysis, we discuss 
a curious behavior of the transverse matrix element in the vicinity of crossings specific to the Hubbard 
Hamiltonian. Finally, to illustrate the connection between level crossings and integrability we consider 
in some detail a simple case of integrable 3x3 matrix systems (Section 6). In this case we conclude 
that there are no "accidental degeneracies", i.e. all degeneracies may be associated with dynamical 
conservation laws. 

The majority of the results do not depend on the number of sites {N) or electrons (M). However, 
cases of even and odd N and M have to be treated differently, the latter case being substantially sim- 
pler because of less symmetry that needs to be taken into account. We therefore take for concreteness 
that both these numbers are even, N = 2n and M = 2m with an equal number of spin up and spin 
down electrons. Throughout the paper analytical results are illustrated by numerical computations 
for a special case of benzene (m = n = 3), and a few other choices of (m, n). 



2 U— independent Symmetries 

The starting point of our analysis of the spectral properties of the Hubbard Hamiltonian ( |1 . 1|) is the 



list of its U-independent symmetries. We will follow closely the paper by Heilmann and Lieb ITT 
adopting in most cases their notation. 

The role of the U-independent symmetries is twofold. First, they greatly simplify the diagonal- 
ization of the Hamiltonian and integrals of motion and provide a convenient labeling for various parts 
of the spectrum. Besides, for the purposes of studying specifically the effect of parameter dependent 
conservation laws all these symmetries have to be factored out anyway. 

U-independent symmetries fall into three major categories - the symmetry of the polygon, the spin 
symmetry, and the particle-hole symmetry. Spatial symmetries can be generated by two operators: 
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operator C2n that rotates the polygon by vr/n, and operator a that reflects it in a line through vertices 
n and 2n (see Fig. |^). These operators can be expressed in terms of on site creation and annihilation 
operators as 

n— 1 n 
C2n = ^ = n n Jfc«;-A:s ^' = II 11 Jfe-l,s;2n-fc,s (2.2) 

s=T,ifc=i s=T,ife=i 
where Jjsj's' interchanges the orbitals (js) and {j's') 

Jjs,j's' ~ 1 ft-js fl-j's' CjgCj'g/ -|- c}jigiCjs (2-3) 

Eigenvalues of C*2n are e*'^^/", where P is an integer < P < 2n — 1 that represents the total 
momentum of the state. The reflection operator a has eigenvalues ±1. 

i = i j = 2 

j = 2n = 6/ \ j = n = 3 





j=n+l=4 



Figure 2: Regular hexagon {n = 3). Cq rotates by 7r/3. a reflects in a line through vertices 3 and 6. 
The spin symmetry can be generated by and 

n — fi 

S^ = ^^ ^+ = (^-r = EsV.i (2-4) 
= ^-^+ + ^+^- + 42 (2.5) 

where 

2n 

ns = J2^js (2-6) 

i=i 

To describe the particle- hole symmetry we first define operators J^°^ and J^^'' ■ The operator J^°^ 
changes the sign of the wave-function each time there is a spin s electron on an odd site, while the 
operator J^'^^ interchanges holes and particles for a spin direction s. 

Ji"^ = n (1 - 2^2,+!,.) ) = n(c], + c,,) (2.7) 

j=0 3=1 

It is convenient to introduce the following combinations of these operators: 

j{o) ^j[o)jio) jW^jWjih) jio) ^jio)jih) z^ = jio)jW ^2.8) 



Note the action of various symmetries on the creation and annihilation operators 

Z\c,^Z^ = (-l)^c,t Z\c,,^Z^ = -cl (2.9) 
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It follows from equations ( p.9|) and ( p.lOp that 
and 

/'•°-'|m|, m^) = \2n — m-f, 2n — m^) 



(2.10) 
(2.11) 
(2.12) 

(2.13) 
(2.14) 



where |m|, m^) is a state with spin up electrons and spin down electrons, and 2n is the number 
of sites. Thus, operators 1^°^ and Z-^ conserve the number of particles only at a half filling with an 
equal number of up and down electrons. 

Equations ( |2.9| ) and ( |2.1CI| ) are used to evaluate various commutation relations for 1*^°^ and Z^. 



{Z|,if} = [a,i/]=0 

i^°^C2n - C2J(°)(-1)" = Z^C^n + C^nH-lf' = Z{a - (-1 



\n-l 



crZt = 



(2.15) 

(2.16) 
(2.17) 



where {A, B} = AB + BA, n = + n^, and T is the kinetic energy operator 



N 



T — T^^ (c^^Cj+i s + c^j+i gC 

j=i s=n 



■jsj 



Z| is a unitary operator that anticommutes with the Hamiltonian. Therefore it defines a mapping 
from the set of operators that commute with the Hamiltonian onto itself, i.e. [A, H] = implies 
[Z^AZ-^, H] = 0. Using this transformation, we can construct a new su{2) algebra (ry-pairing su{2)) 
from the spin su{2) 



12, 13, 14, 1 



nt + ni 



2n 



— n I/. 



The Casimir operator of this new 5^(2) 



z^ s'^ Z^ 



(2.18) 



(2.19) 



preserves the number of particles and therefore its eigenvalues L{L + 1) can be included in the set 
of quantum numbers used to label the eigenstates. The combination SU{2) x SU{2)/Z2 yields the 
complete 5*0(4) symmetry of the Hubbard Hamiltonian jl^ , |13|, |14|, p!5| . 
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3 Conserved Currents 



Now let us turn to the parameter dependent integrals of motion. In principle, an infinite number of 



these integrals can be obtained by methods outlined in |T^, |T7|, 0. Clearly, on a finite lattice only a 
finite number of these integrals are independent (see Section 6 for more details). The general form of 
the rth conserved current is 



fc=0 



where are parameter independent operators. The U-independent part of Ir has a simple form 

N 

liU = 0) = = ±{i)^+^T' E (cksc,. - {-irc],c,^r,s) (3.20) 

As far as permanent degeneracies are concerned, an important feature of the currents is that odd 
(r = 2/c + 1) and even (r = 2k) currents transform differently (equations (p.23|) -( ^.26|) ) under spatial 
reflections (a) and partial particle-hole transformation {Z^). Odd currents are in many respects 
similar to the Hamiltonian, while properties of even currents are essentially different. 



First few nontrivial currents were derived explicitly in [0, |T9|, ( See also |^ and ||2^ for 

a discussion on the derivation and the structure of the higher conserved currents.) To analyze the 
benzene example we will need the explicit form of only the first two nontrivial integrals of motion. 
Since Ii = H, these two are 



N N 



h = -iTj2J2 (Ci+2sCi. - c],Cj+2s) -^UJ2Y. (S+lsCj. - c},Cj+i,) (%+!,_, + - 1) (3.21) 
i=i s=n i=i s=n 



N N 



i=i s=n i=i s=u 

N 



j 1 •^~TJ- 



^ ^ 1\ 1 



nn - 2; [^n - 2 



(3.22) 

Odd currents have the same U-independent symmetry as the Hubbard Hamiltonian, i.e. in all 
commutation relations of Section 2 the Hamiltonian can be replaced with any other odd current. For 
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even currents two of the commutation relations are different. Namely, for odd currents one can show 
that 

['T,Wi] = (3.23) 

{Z^,/2fc+i} = (3.24) 

while for even currents 

{a,/24 = (3.25) 

[Z^,J2d = (3.26) 

In other words, equations ( p.23| - pl2^ ) mean that even and odd currents have different a and Z| 
parities. For k = 1 equations ( p.23| - pl26D have been derived in [^. Note also that since both a and 
Z| are unitary and Ir is Hermitian, equations ( p.25| ) and ( 3.24 ) imply 



Trir = 

Another property of conserved currents which we will use bellow is 

j(°)/,(T,f/)j(°) = {-lY^'iri-T,U) 



(3.27) 



(3.28) 



Equations ( |3.23| )- (|3.26| ) and (|3.28| ) for I2 and Is can be verified by inspection. For higher currents 
they can be derived by establishing the transformation properties of the transfer matrix obtained in 

Finally, let us measure U, T and all other energies in units of f/ — AT. This is equivalent to the 
replacement 

U-AT=1 U = u T = {u~ l)/4 (3.29) 

We see from ( P . 28|) that the spectrum and the eigenfunctions of Ir for T < and T > are related 
via a simple transformation and it is therefore sufficient to consider T < 0. Thus, we have to compute 
the spectra only for < u < 1. 



4 Diagonalization 

In this section we outline a method that we use for detailed numerical study of the spectra of the 
Hamiltonian and conserved currents. In subsequent sections we will compare the computer generated 
results for benzene (m = n = 3) and several other values of m and n to general predictions based on 
symmetry. 

From the results of Section 2 it follows that the set of quantum numbers that label the eigenstates 
of the Hamiltonian and integrals of motion is {P, a, S, L}. However, there are certain restrictions 
on this set since not all symmetries and integrals of motion mutually commute. The quantum number 
1^°^ can be assigned only when m = n and a can be used only when P = or P = n (see ( 2.14 ) and 



( [4.34D ). Also, since a changes P to 2n — P without affecting the energy and other quantum numbers. 



the spectrum for 2n > P > n is the exact copy of that for n > P > 0. Therefore, it is sufficient to 
diagonalize the Hamiltonian and the currents for P ranging from to n. Finally, since even currents 
anticommute with a, their eigenvalues cannot be specified whenever a is specified. 
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To utilize the conservation of momentum, we introduce in the usual way the momentum space 
creation and annihilation operators and dps 



-| 2ra 1 2n 

< = ^ E ^""''^'^ dps = ^Y: e-^'^^^/"c,. (4.30) 
V2n =1 v2n ,=i 



where p = 0,1, . . . ,2n — 1. 

Next, we choose a basis for 2m-electron wave functions with m spin up and m spin down electrons. 

|r) = |p;q)=<t...4^T4a---4.JO) (4-31) 
with the following ordering convention 



Pi <P2< ■■■ <Pm 
qi<q2< ... <qm 



(4.32) 



Components of |r) that are not in the interval [0, 2n — 1] have to be reduced modulo 2n into this 
interval. If the components of a vector are not in the order ( [4.32|) they have to be permuted to obtain 
these ordering and the wave function should be multiplied by (—1) if the permutation is odd. The 
total momentum of the state is 

2m 

P = ^ Tfc (modulo 2n) (4.33) 

k=l 

From ( p.9| - p.H"D and (^4.30|) we identify the action of operators of U-independent symmetry on 
basic states (|4.31| ) 

(T|r) = |-r) (4.34) 

jW|r) = (-l)2™+^'|_(r + ne)) (4.35) 

Z^|p;q) = (-l)(™+«)|(p + ne);-q) (4.36) 

where Q = Z^feLi Qk is the total momentum of spin down electrons, r denotes taking the complement 
of the set TZ = {vk} in the set of integers from to 2n — 1. Finally, +ne stands for adding n to each 
component of a vector. 

The rotation operator commutes with the Hamiltonian, conserved currents, S^, and L^. 
Therefore, these operators split into blocks corresponding to different values of P. This significantly 
reduces the size of the matrices to be diagonalized. To generate the blocks for each value of P we 
need to derive the matrix elements of the Hamiltonian, integrals of motion, and symmetry operators 
in the basis ( |4.31| ). Equations ( ^.34| ) and ( ^.35| ) allow to write down a, and Z-^ in the matrix 



notation. Matrix elements of the Hamiltonian, and the operators S^, I2, and I3 are summarized in 
Appendix B. 

Now we have all ingredients needed to produce a computer program (e.g. on Mathematica) that 
generates exact blocks of the Hamiltonian, symmetry operators, and first two currents for any value 
of the total momentum P. 

After the value of the total momentum P is chosen, the program picks one of the remaining 
U-independent symmetries, say a, and finds a unitary transformation that makes it diagonal. The 
columns of the matrix of this transformation are the eigenvectors of a. Since we know the eigenvalues 
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of 6" exactly (±1), determining the eigenvectors for each eigenvalue reduces to a set of linear equations 
which can be solved analytically. 

In this way, step by step, we diagonalize all parameter independent symmetries and split the 
Hamiltonian and conserved currents into smaller blocks. If the system is not too large, this can 
be done analytically using Mathematica. For example, the complete Hamiltonian for benzene is a 
400 X 400 matrix, while individual blocks range in size from 1 x 1 to 16 x 16. 

Finally, when all U-independent symmetry is exhausted, we get exact blocks of the Hamiltonian 
and conserved currents with all U-independent symmetry quantum numbers assigned to each block. 
Each block is then diagonalized individually at different values of u ranging from to 1 (see ( t3.29|) ). 
Often, the small size of a block allows for analytical diagonalization. 

5 Permanent Degeneracies and Other Spectral Properties 

As soon as the Hamiltonian (|1 . 1| ) is diagonalized, we discover numerous cases of permanent degeneracy, 
see e.g. Table |l] and Fig. [|, |^, and ^ We make the following observations. The majority of degeneracies 
are twofold degeneracies with respect to cr = ±1. However, there are also several multiply degenerate 
levels. The characteristic feature of these states is that the energies are linear in U. 



p 


n 


m 


(degree of degeneracy, number of levels) 





3 


2 


(1, 12) 


(2, 12) 


3 


3 


2 


(1, 15) 


(2, 7); (3, 1); (7, 1) 





3 


3 


(1, 32) 


(2, 18) 


3 


3 


3 


(1, 30) 


(2, 12) 


(7, 2) 





4 


2 


(1, 28) 


(2, 31) 


(3, 1); (7, 1) 


4 


4 


2 


(1, 15) 


(2, 41) 


(3, 1) 





4 


3 


(1, 70) 


(2, 150); (3, 2); (8, 2) 


4 


4 


3 


(1, 56) 


(2, 155); (3, 1); (7, 1); (8, 2) 





4 


4 


(1, 126); (2, 229); (3, 1), (7, 1); (12, 2) 


4 


4 


4 


(1, 126); (2, 229); (3, 1); (7, 1); (12, 2) 



Table 1: Degrees of degeneracy in the sectors P = {0,n} for n = {3,4} 



First, we analyze the twofold degeneracies. Let ipi^E, hk+i, cr, S, L), where E is the energy, be an 
eigenstate of the Hamiltonian as well as all other odd currents and operators a, S^, and L^. Being an 
eigenstate of a, the state ip cannot be an eigenstate of any even current since l2k and a anticommute. 
Let 4>{E, l2k, P, S, L, l2k+i) denote an eigenstate of all even and odd currents. On the other hand 
(p = ip when a is not assigned, i.e. when P ^ {0, n}. We note that equation (|3.25|) means that either 
the state ip is annihilated by the operator l2k, l2k'ip = for all U, or 



l2ki'{E, hk+i, cr, S, L) = tp{E, hk+i, -cr, S, L) 
Equations ( p. 26] ) and ( |3.25| ) imply 

Z|V^„(E, /2fc+i, S, L) = ipni-E, -hk+i, L, S) 



(5.37) 



(5.38) 
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Zi(j)n{E, hk, S, L, hk+i) = <j)n{—E, l2k, L, S, —hk+i) (5.39) 
aME, hk, P. S, L, hk+i) = ME, -hk, -P, L, hk+i) (5.40) 

It follows from equation (|5.38|) that for each state of energy E and eigenvalues of odd currents {l2k+i} 
there is a state with the energy —E and eigenvalues {—I-zk+i}- According to equation (|2.13|) these 



two states have the same number of particles only ii m = n. Equation ( |5.39|) implies that eigenstates 



of l2k are doubly degenerate. Finally, from ( |5.37| ) we conclude that any state in the sector P = {0, n} 
that is not annihilated by all l2k is at least twofold permanently degenerate. This applies to all values 
of m and n. In other words, all energy levels in the sector P = {0, n} that are not in the kernel of l2k 



for all U and k are doubly degenerate (see also pO[). Consequently, all nondegenerate states should 
be annihilated by any even current. 

A twofold degeneracy is explained on symmetry grounds as soon as we identify an even current 
that maps the two degenerate states into each other. Unfortunately, we do not yet know how to 
prove for arbitrary m and n that for any twofold degeneracy there exists such an even current. For 
benzene we checked numerically that all doubly degenerate states are maped into each other by I2 
and therefore no other even currents are needed to explain twofold degeneracies. 

Since nondegenerate states are annihilated by any even current, nondegenerate states exist only 
if all even currents have nontrivial kernels. Let us show that this necessary condition is met when m 
is odd 01 m ^ n. Indeed, since l2k maps the subspace a = +1 into the subspace a = —1 ( ^.37|) , its 
kernel has to contain at least \d^i — d^i\ states, where d±i are the dimensions of a = ±1 subspaces. 
Since the eigenstates of cr = ±1 are |r) ± |— r) ( [4.34| ), d^i — d^i is the number of states such that 



|r) = |— r) minus the number of states such that |r) = — |— r). This is calculated to be 

^.^-^-^ = i("rJ]H'r^t-J]y (5.41) 



where [x] is the integer part of xQ. 

The evaluation of an upper bound on the fraction of the nondegenerate states (/) in the sector 
P = {0, n} can be reduced to a combinatorial problem. We note that the number of such states 
cannot exceed the dimension of fcer/2 at f/ = restricted to the subspace P = {0,n}. Let us denote 
this dimension by K. The kernel of /2 at [/ = consists of all states |r) such that 



E 

k=l 



. 2TTpk , . 2TTqk 

sm h sm 

n n 



(5.42) 



Thus we have / < K/D = g, where D is the dimension of P = {0,n} subspace. First, one can show 
by explicitly constructing a P = state such that the left hand side of (|5.42| ) is nonzero that K < D 
for 2n — 2 > m > 2 and n 7^ 2. This means that twofold degenerate states exist for all values of m 
and n except n = 2 and m = {0, 1, 2n — 1, 2n}. Numerical values of g for the first few m and n are 
tabulated in Table 0. Since g{m, n) = g{2n — m,n), only values for m < n are shown. 

We see from Table |I] that at small m and n a substantial fraction of levels are nondegenerate. 
Next, we consider the limit n ^ m. The exact value of D depends on the greatest common divider 

^The high degeneracy of the zero eigenvalue of even integrals of motion suggests that even conserved currents have 
an additional symmetry. An interesting open question is whether the corresponding eigenstates can be chosen to be 
independent of U (compare to the discussion on multiply degenerate levels of the Hamiltonian bellow). 
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m\ n 


2 


3 


4 


5 


6 


2 


1.0 


0.89 


0.79 


0.76 


0.71 


3 




0.88 


0.765 


0.70 


0.64 


4 






0.764 


0.675 


0.62 


5 








0.670 


0.610 


6 










0.608 



Table 2: Values of K/D. Note that the upper bound on the number of nondegenerate states mono- 
tonically decreases as m and n increase. 

of m and n. However, this is not essential for large n, since we can always decrease or increase n by 
a small number so that m and n become mutually prime, in which case the dimension of P = {0, n} 
subspace is 1/n oi the total dimension of the Hilbert space. 

,„2m-l 




^2m-i for 72 > m (5.43) 

To analyze the behavior of K in the limit n ^ m, we note that there are two ways how a number of 
terms in ( ^.42 ) can cancel. The first option is the pairwise cancellation. Alternatively, r terms can 



cancel if their momenta correspond to the vertices of the regular r-gon. This can happen only if n is 
divisible by r. For a given p there are at most four other momenta with which it can cancel pairwise. 
Therefore, for n ^ m the number of ways in which the sum ( |5.42| ) can cancel pairwise grows slower 



than n™, while the number of non-pairwise cancellations grows slower than n*" ^. We conclude that 

„ ^ const „ ,r AA\ 

f < for n > m (5.44) 

This asymptotics combined with numerical values from Table ^ suggests that almost all states in the 
sector P = {0,n} are at least twofold degenerate in the thermodynamical limit. 

As for the multiply degenerate states the first observation is that, at least for benzene, all of 
their eigenvalues are linear in the coupling parameter u. This seems to suggest that these states are 
simultaneous eigenstates of the kinetic (T) and the potential (U) energy operators (|3.29|) . Indeed, this 
turns out to be the case for the majority of the multiply degenerate levels for benzene. However, the 
benzene Hamiltonian also has two eigenstates with E = ±u/2 and P = 3 that are not simultaneous 
eigenstates of T and U. 

Some of the multiply degenerate states can be obtained using the prescription of |l^. The idea 
is the following. We take a state |r) that has only one species of electrons, say 21 (/ < n) spin down 
electrons, and no spin up electrons 

|r) = |0;q) = |0;gi,...g2i) (5.45) 
Clearly, |r) is an eigenstate of both T and U. 

f|r) = ^^f]cos^|r) U\r) = u{n/2 - l)\r) (5.46) 
2 k=i ^ 
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On the other hand, |r) is a lowest weight state in the so(4) algebra constructed from S and L 

S,\r) = -l\r) S = l 



(5.47) 

4|r) = (/-n)|r) L = n-l 



The energy according to (|5.46|) is 



= ^ cos — + M(n/2 - /) (5.48) 

2 k=i ''^ 

If m > /, we can get an eigenstate |f) with m spin up and m spin down electrons, by applying the 
raising operators of the so(4) algebra 

I?) = ^^L™-'|r) (5.49) 
The total number of states that can be obtained in this way is 

M^gg) (5,50, 

f2n\ 

while the number of states that have the same energy ( |5.48|) cannot exceed r,; • For m = n = 3 the 



allowed choices for I are / = {0, 1, 2, 3}. The choice / = 1 yields 15 eigenstates of which 3 states have 
the total momentum P = 3 and are degenerate with the energy E = u/2. By partial particle- hole 
symmetry (|5.38|) the choice / = 2 also gives 3 states with P = 3 and the energy E = —u/2. The 
actual number of states of energy u/2 and P = 3 is 7 (Table |l|). Wave functions of 6 of these states 
are independent of U. 

We arrive at the conclusion that contrary to what was conjectured in |l^, not all U-independent 



eigenstates can be obtained using the above prescription. The six U-independent eigenstates together 
with the corresponding states that have lowest weight in so(4) are given in Appendix A. 

Clearly, the multiple degeneracies discussed above mean that there is an additional U-independent 
symmetry of the Hubbard Hamiltonian that have been overlooked in Section 2. This symmetry seems 
to be related to the high degeneracy of the spectra of the potential and kinetic energy operators. 
However, currently we do not know what is the algebra of generators of this additional symmetry and 
how its representations are to be classified. 



12 



0.75 ■ 



-0 . 75 ■ 



-Q . 6 



0.8 



Figure 3: Energies for (P, a, /(°), ^, L) = (0, 1, 1, 1, 1) (cf. Fig D in units of f/ - 4T as functions of 
u = U/{U — AT). Levels marked with crosses are twofold degenerate with respect to a = ±1. They 
also correspond to {P, a, S, L) = (0,-1,1,1,1) (equation (|5.37|) ). Note that since S = L the 
spectrum is symmetric under the reflection E —>■ —E (equation (|5.38|)). 



\l3 



0.05 



■■:-6. . 



Figure 4: Eigenvalues of Jg for (P, a, /(°\ 5, L) = (0,1,1,1,1) (cf. Fig. |) in units of ?7 - 4T as 
functions of u = U/{U — 4T). All eigenvalues are twofold degenerate with respect to a = ±1. They 
also correspond to {P, a, S, L) = (0,-1,1,1,1) (equation ( |5.37| )). The spectrum is symmetric 
under the reflection Jg — Jg (equation (|5.38|) ). Note however that crossings of J3 and H occur at 
different values of u except u = {0, 1}. 
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u 



-0.5 ■ 



Figure 5: Energies for (P, a, /(°), 5, L) = (0, -1, -1, 0, 1) (cf. Fig. | and 0) in units of U - AT as 
functions oi u = U/{U — AT). Levels marked with crosses are twofold degenerate with respect to 
a = ±1. They also correspond to {P, a, S, L) = (0,1,-1,0,1) (equation (|5.37| )). Upside down 
the figure shows levels for the same set of quantum numbers with S ^ L (equation ( p.38| )). Note the 
nondegenerate levels. 
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-0.2 



-0.3 



■■0-,.,2 



4 ::g:;::fi: 



Figure 6: Eigenvalues of Is for (P, a, &\S, L) = (0, -1, -1, 0, 1) (cf. Fig. | and 0) in units of U- AT 
as functions of u = U/{U — AT). Levels marked with crosses are twofold degenerate with respect to 
a = ±1. They also correspond to {P,a, S, L) = (0,1,-1,0,1) (equation ( |5.37| )). Upside down 
the figure shows levels for the same set of quantum numbers with S ^ L (equation (|5.38|) ). 
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Figure 7: Eigenvalues of J2 for (P, S, L) = (0, —1, 0, 1) in units oi U — AT as functions of m = 
U/{U — AT). Note the difference as compared to Fig. ^ and ^ The quantum number a cannot be 
assigned since I2 anticommutes with a ( p. 25 ). All eigenvalues are twofold degenerate with respect to 
S ^ L and E —E (equation (|5.39| )). Since P = —P = the spectrum is symmetric under the 
reflection (equation ( p.40| )). 




-1.5 



Figure 8: Energies for (P, S, L) = (2, 1, 0, 0) (cf. FigPandp!0|) in units of U—AT as functions of m = 
U/{U — AT). Since S = L the spectrum is symmetric under the reflection E —>■ —E (equation ( |5.38|) ). 
Note the level crossings. 
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-1 



Figure 9: Eigenvalues of J2 for {P, S, L) = (2,1,0,0) in units of U — AT as functions oi u = 
U/{U — 4T). All eigenvalues are twofold degenerate with respect to E ^ —E (equation p.39D . Note 
that unlike Fig ^ and |10| there is no reflection symmetry since P 7^ —P (equation (|5.40| )). 



0.1 



-0.1 



Figure 10: Eigenvalues of J3 for (P, &\ S, L) = (2, 1, 1, 1) (cf. Fig | and D in units of U - AT as 
functions of u = U/{U — AT). Since S = L the spectrum is symmetric under the reflection E —E 
(equation (|5.38|)). 
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6 Level Crossings 



The spectra of the Hamihonian and conserved currents display numerous cases of crossings of levels 
that have the same U-independent symmetry quantum numbers. In this section we argue that this 
type of behavior is in fact expected in integrable models. 

Let us first clarify the notion of a nontrivial integral of motion for a quantum system on a finite 
lattice. Let H{u) = Hq + uV be a Hamiltonian that, in a certain basis, can be represented by an 
s X s matrix. For example, H{u) can be one of the blocks of the Hubbard Hamiltonian. Clearly, 
any matrix M that is an analytical function of H{u) and u, M = f{u,H{u)), commutes with H. 
Obviously, M is not an independent conserved quantity. Therefore, we say that J{u) is a nontrivial 
integral of motion if J{u) is hermitian, [J{u), H{u)] = and at the same time J{u) cannot be written 
as an analytical function of H{u) and u. 

Now let us show that the existence of a nontrivial integral of motion implies at least one level 
crossing. Since J{u) and H{u) commute, there is a basis where both these operators are diagonal. 
Let Ei{u), . . . , Es{u) and Ji(u), . . . , Js{u) be the eigenvalues of H{u) and J{u) respectively. Consider 
the following set of algebraic equations: 

aiEl~^{u) + . . . + as-iEi{u) + as = Ji{u) 



(6.51) 



aiEl '^{u) + ... + as-iEi(u)+as = Js{u) 



If equations ( |6.51| ) have a solution, J{u) and H{u) are not independent. Namely, J{u) = aiH^~^ + 
. . . + tts- The system (|6.51|) has no solutions if and only if Ei{u*) = Ej{u*) and Jiiu*) ^ Jjiu*) for 
some i, j, and u* . Thus, the Hamiltonian can have a nontrivial integral of motion only if it has a 
level crossing. 

Let us analyze the suppression of level repulsion in integrable models in more detail. Here we 
consider only real Hamiltonians that depend linearly on the coupling constant H{u) = Hq + uV . We 
assume that there is a value of u = tt at which all eigenvalues of H{u) are nondegenerate. This is 
true for instance for blocks of the Hubbard Hamiltonian for benzene. In general, since all integrals 
of motion mutually commute, permanent degeneracies occur only with respect to quantum numbers 
of parameter independent symmetries. Therefore, after all these symmetries are factored out, at a 
certain value of u all states in a given block are nondegenerate. Hence, we can use the nondegenerate 
perturbation theory to write down the variation of the eigenvalues and the matrix elements of the 
perturbation with u in the vicinity of u. 



du 



Vnn{u) (6.52) 



dV^ = 2Y.^r^P^, (6.53) 

du f^^En{u)-E,{u) 

du En{u) - Ej{u) .^^ Era{u) - Ej{u) 



where all matrix elements are evaluated in the running basis 

Htpn{u) = En{u)ll)n{u) 

Vi,{u) ^ {Uu)\V\ij^{u)) 
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Equations ( |6.52| )-( |6.54| ) were used in to derive the distribution of energy eigenvalues of the 
irregular spectrum in the semiclassical limit. Note in particular that these equations are not model 
specific. The model itself is defined by the initial conditions - the values of Vij{u) and Ei{u). If 
the initial conditions in any block are slightly perturbed^, the integrability is lost and the majority 
of crossings are converted into anticrossings (see e.g. Fig. [TI| and At the same time, since the 
perturbation has the block diagonal structure of the original Hamiltonian, the perturbed Hamiltonian 
still has the same U-independent symmetries. 

Now consider a pair of levels that get close at m = Mq, i.e. we assume that the absolute value of 
the energy difference between these two levels A(u) = Ei{u) — E2{u) is much smaller than the energy 
distance to the remaining levels. Then, from (|6.52|) and ( |6.53|) we obtain the following equations for 
A(u): 



dA 

du 
d^A 

du 



22 



A 



(6.55) 
(6.56) 



where 



F(u) 



2E 



EAu)-EAu) 



E, 



u 



We can interpret Ei and E2 as coordinates of two one dimensional particles. Particles move with 
a relative velocity (Vii — V22) and interact with a force 4yj^|/A + F{u). If V12 7^ as A —> 0, an 
infinite repulsion prevents the particles from colliding. On the other hand, if levels 1 and 2 have 
different U-independent symmetry, Vu^u) is identically zero and levels are permitted to cross. These 
arguments constitute the essence of the noncrossing rule. 

Let us incorporate conservation laws into the picture. Let the corresponding block of the conserved 
current be J{u*) = J{u) + {u* - u)W + (m* - u^W' + . . . Since J{u*) commutes with H{u*) to all 
orders in {u — u*), we get 



[H{u*), J{u*)] = ^ [H{u), W] + [V, J{u)] = 
Evaluating the matrix element of the above equation between states 1 and 2, we find 

E^{u)-E2{u) 



(6.57) 



VM = Wu{u) 



Jliu) - J2{u) 



const A — s> as A 



(6.58) 

and levels are 



Therefore, unless Ji{u) = J2{u) in the vicinity of -u = uq, V12 
allowed to cross. 

We have checked for benzene that the eigenvalues of the Hamiltonian, I2, and I3 are never degen- 
erate at the same value of u except at two points of special symmetry u = and u = 1. Because I2 
also connects all doubly degenerate states, we conclude that I2 and I3 together with U-independent 
symmetries explain all cases of level crossings and permanent degeneracies in the spectrum of the 
Hubbard Hamiltonian for benzene. 



^e.g. by adding small random numbers to Ei{u) and Vij{u). 
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-1.5 



Figure 11: The block for {P, S, L) = (1,1,0,1) has been perturbed by adding small random 
numbers to Vij{u) and Ei{u) at tt = 0.6. Energies at other values of u are determined according to 
( |6.52[ )-( p.54| ) and plotted in units of U — AT as functions of u = U/{U — 4T). 



1.5 



0.5 



-0.5 



-2 L 



u 



0,2 0--4::=;;;!=».»=::::...Q„„.6,.,. . 



Figure 12: Unperturbed energies for (P, S, L) 
u = U/{U -AT). 



(1, 1, 0, 1) in units of U — AT as functions of 
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In the vicinity of a crossing at u = u* the transverse matrix element Vi2{u) can be expanded in 
series in the energy difference A{u). 



yi2(A) = aiA + asA^ + . . . 



(6.59) 



The coefficient at the hnear term can be derived using the degenerate perturbation theory. 

a, = y, (6.60) 



Ei{u*) 



Here the basis at m = m* has to be chosen so that Vuiu*) = (see e.g. Chapter 4 in 

The surprising feature of the Hubbard Hamiltonian is that, at least for benzene, the linear term 

vanishes for all crossings at m 7^ {0, 1}, while individual terms in the summation ( |6.60| ) do not. At 

this point, we do not have an explanation for this seemingly puzzling phenomenon. 

To see that this is not a mere consequence of integrability, we study a simple example of 3 x 3 

Hamiltonians. This example also illustrates the connection between crossings and integrability and 

is therefore interesting on its own right. 
Let 

H = Ho + uV I = Io + uW (6.61) 
be two 3x3 real symmetric matrices. We call a pair {H, I) integrable if these two matrices commute 



[H,I] = 



(6.62) 



Since we can always add multiples of the identity matrix to H and / without affecting integrability 
( I6.62D or any of the level crossings, with no loss of generality, we can assume that both H and / are 
traceless 

TrH = TrI = 

Further, we say that a pair {H, I) is trivial if the u dependence can be eliminated from either H or I 
by a u-independent unitary transformation (change of basis) and (or) by taking linear combinations 
of H and /. 

First we demonstrate that if {H, I) is a nontrivial integrable pair, both H{u) and I{u) have a 
single level crossing. Indeed, one can check that by changing the basis and taking linear combinations 
any nontrivial integrable pair can be brought to the following "canonical" form: 



(6.63) 





" 1 















Vl3 


H = 





1 





+ u 





V22 


V23 










-2 






V23 


-Vn - V22 _ 




' 1 





" 




Wu 


Wl2 


Wis 


I = 





-1 





+ u 




W22 


W23 















Wi3 


W23 


-Wn - W22 



(6.64) 



Matrix elements of W can be written in terms of fj/s using ( |6.62|) . 

5vl 



Wn 



f iif 22 + 2wf3 



4^22 



4^23 



9{vn - V22) 



(6.65) 
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wi2 = ^'"^^^^^ ^y^3 = ^^3/3 (6.66) 
3{vii - V22) 

9{vn - V22) 

W23 = -W23/3 W33 = -wu - W22 (6.68) 
Since W has five independent matrix elements, while V has only four, Wi/s are constrained by 

W33 + ^ h T -W13W23 = (6.69) 

3W23 3wi3 3 

-?/(«) has an apparent crossing at m = 0. Further, it turns out that the matrix I{u) has a crossing 
if and only if the constraint ( |6.69D is met! This can be verified either directly or using where. 



among other things, degeneracies of real symmetric 3x3 matrices are analyzed in detail. The crossing 
of / is at 

2/u* = {W22 - Wu) + Wi2iwi3/W23 - W23/W13) (6.70) 

Similarly, it can be shown that linear combinations of H and / always have crossings too. Therefore, 
if {H, I) is a nontrivial integrable pair, both H{u) and I{u) have a single level crossing. 

Now let us show that the converse is also true, i.e. that any real symmetric 3x3 matrix that has 
a pairwise crossing is also integrable. We start with a real symmetric 3x3 matrix H = Hq + uV that 
has a pairwise crossing aX u = u* . By redefining the parameter u ^ u — u* and a suitable choice of 
basis we can write H in the form ( |6.63|) . Therefore, there exists a matrix I of the form (|6.64|) such 



that [I-,H\ = 0. In view of this explicit construction of a "dynamical conservation law", we may say 
that there are no accidental degeneracies in the case of 3 x 3 real matrices. 

Finally, we note that the linear term in expansion (|6.59| ) is non zero for both H and / unless 
specially arranged. 



7 Conclusion 
7.1 Summary 

We have shown that for all m and n except m = {0,1, 2?t, — l,2n} and n = 2 there are twofold 
permanent degeneracies for the values of the total momentum P = and P = n. These degeneracies 
are a consequence of different transformation properties of even and odd currents with respect to 



spatial reflections (a) and partial particle- hole transformation (Z|) (equations (|3.23| - p^26D ). We have 
argued that in the thermodynamical limit the fraction of doubly degenerate states in the sector 
P = {0,n} approaches one. For benzene (m = n = 3) we have checked that all doubly degenerate 
states are maped into each other by the first even current l2- We have also seen that some states can 
be nondegenerate since all even currents have nontrivial kernels. 

We have learned from the benzene example that there is a number of multiply degenerate states 
that are simultaneous eigenstates of the potential and kinetic energy operators. The wave functions 
of these states do not depend on the coupling U. This suggests the existence of an additional U- 
independent symmetry. However, as we have seen in Section 5 and Appendix A some of these states 
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cannot be obtained by applying the raising (lowering) operators of the so(4) algebra to states that 
have only one species of electrons. 

Further, we have demonstrated how the constraints of the noncrossing rule are lifted in integrable 
models. In fact, we have argued that in integrable models crossings of levels with the same parameter 
independent symmetry quantum numbers are expected. This was also illustrated by a simple example 
of 3 X 3 matrix models where there is a one to one correspondence between crossings and integrability. 
We have seen that for m = n = 3 the first two nontrivial conserved currents are sufficient to explain 
all cases of level crossings. We have shown numerically that in the case of benzene the transverse 
matrix element is quadratic in energy difference in the vicinity of all pairwise level crossings. 



7.2 Open Questions 

1. What is the symmetry responsible for the multiply (more than twofold) degenerate levels in 
the spectrum of the Hubbard Hamiltonian? Are energies of multiply degenerate states always 
linear in the coupling parameter? How one can count the degrees of degeneracy and the number 
of multiply degenerate states for arbitrary m and n? 

2. What is the symmetry responsible for the high degeneracy of the zero eigenvalue of even integrals 
of motion? 

3. For given m and n what is the exact number of states such that /2fc^ = for all k and Ul 
What is the exact number of twofold degenerate states with P = {0, n}? 

4. Is the transverse matrix element always quadratic in energy difference in the vicinity of all 
pairwise level crossings in the Hubbard model? What is the explanation of this phenomenon? 

8 Appendix A 

Here we write down U-independent eigenstates that have E — u/2. See Section 4 for notations. 

|ri) = -|0,1,2;1,2,3)-1/2|0,1,3;0,2,3) + 1/2|0,3,4;0,3,5)+ 

|0, 4, 5; 3, 4, 5) - 1/2|1, 2, 4; 1, 2, 5) + 1/2|1, 4, 5; 2, 4, 5) + (T^i) 



S 



L 
-0 U 



2 

: 



(8.71) 



|r2) = |0,2,5;1,3,4) + |0,3,4;1,2,5) + |0,3,5;1,2,4) 

-1 1, 4, 5; 0, 2, 3) - |2, 3, 5; 0, 1, 4) - |2, 4, 5; 0, 1, 3) + (T^i) 



S^O L = 1 



(8.72) 
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Irs) = -|0,2,5;1,3,4) + 1/2|0,3,4;0,3,5) + 1/2|1,2,4;1,2,5) 

1/2|1, 4, 5; 2, 4, 5) + 1/2|0, 1, 3; 0, 2, 3) + |0, 1, 4; 2, 3, 5) - (T^i) ^ ^ = ° 



|r4) = -|0,1,3;0,2,3)-|0,3,4;0,3,5)+ 

|l,2,4;l,2,5) + |l,4,5;2,4,5) + (t^i) = ^ = ° 



Irs) = |0,1,2;1,2,3) + |0,1,4;2,3,5)|0,1,5;2,3,4)+ 

|0, 2, 4; 1, 3, 5) + |0, 2, 5; 1, 3, 4) + |0, 4, 5; 3, 4, 5) - (T^i) " ° " ° 



(8.73) 



(8.74) 



(8.75) 



Ire) = 1/2|0,1,3;0,2,3)+|0,1,5;2,3,4)-|0,2,4;1,3,5)- 

1/2|0, 3, 4; 0, 3, 5) - 1/2|1, 2, 4; 1, 2, 5) + 1/2|1, 4, 5; 2, 4, 5) - ^ ° ^ ° 

(8.76) 

The corresponding lowest weight eigenstates of the Hamiltonian are 

|fi) = -|l;2)-|2;l) + |4;5) + |5;4) = Q = -2 (8.77) 

|f2) = |0,1;2,3)-|0,2;1,3)-2|0,3;1,2)-2|0,3;4,5)-|0,4;3,5) + |0,5;3,4) ^ ^ 

-2|1, 2; 0, 3) - |1, 3; 0, 2) + |2, 3; 0, 1) + |3, 4; 0, 5) - |3, 5; 0, 4) - 2|4, 5; 0, 3) " ~^ 

(8.78) 

|f3) = |;l,2) + |;4,5) ^, = -1 L, = -2 (8.79) 
|f4) = -2|0;3)-2|3;0) + |l;2) + |2;l) + |4;5) + |5;4) S,^0 = -2 (8.80) 
|f5) = |;0,3) S,^-l L, = -2 (8.81) 

|f6) = |;l,2)-|;4,5) ^ -1 = -2 (8.82) 

States Ifs), Ifs), and Ife) describe one species of electrons (spin down). The remaining states 
cannot be obtained by the prescription detailed in the end of Section 5. 
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9 Appendix B 



In this section we collect the expressions for the matrix elements of the Hamiltonian, 12, I3 and S^. 
Let pi < . . . < pt be the momenta of spin up electrons that are in P = {pk} but not in V' — {p^}, 
Pt+i < ... < p2t be the momenta that are in V' but not in V. Similarly, we can define Q — {qk}, 
qi < . . . < Qu, and qu+i < . . . < q2u- The number of elements in the intersection of V with V' and Q 
with Q' is therefore L{V n V') — m — t and L{Qn Q!) — m — u. We also define 

e(r', r) = (r'l^+^T^PiT ■ ■ ■ 42t^dp,^4u+ildqil ■ ■ ■ 42U^9a|r) = ±1 (9.83) 
In this notation the matrix elements of S^, H, I2 and I3 take the following form: 



(r'|S» 



n~L{Vn Q) if |r) = |r') 

e(r', r) if 7^ = 7^' and = i = 1 

otherwise 



(9.84) 



W\H\r) 



(r'|/3|r) 



(r1/2|r) 

(u-l) 



2m 



U-l ^ -KTk. U , 

—— 2^ cos 5r',r - :3-e(r , Y)6uiOnbpp> 



V sin — ^(5r',r - 77-e(r', r)5ui5nbpp' 
4 n 2n 



k=l 

3 2m „,2 



128 

u{u — 1) 



ESvrrfc u (u 1) . -^—^ nj^i^ ^ i^^^ 
^'"^ ^2 — cos 



m Tr{qj - qi) 



k=l 

2 m 



n 



y ^ sm cos Or',r + e(r , r 



k,j,l=l 



n 



n 



(9.85) 
(9.86) 



8 



k=l 



n 



n 



8n2 



k=l 



COS ■ 



n 



cos 



7r(p2 - gi - p'k) 



n 



+ p^q 5y,i5a5pp' + e(r',r)- 



, M^(M-l)r 7r(g4-g2) 



4n2 



(cos ■ 



cos 



7r(g4 - qi) 



n 



Spp' V cos — (5u2(^to + (p -s-^ g, t)]+ 



, m2(m-1) 
7r(g3 - 52 +P2) 



COS + COS + COS ^(g3-g2-pi)^ 



n 



n 



n 



cos 



SuiSt2 + {p ^ q,u ^ t) 



5ppi— 



, u{u-lf 



e(r',r) 



16n 



cos 



7r(pi+p2).„ 7r(gi-g2) 



n 



(2 cos ■ 



7r(pi - qi) 7r(p2 - ^2) 

cos COS ^ — COS 



COS 



n 

7r(p2 - qi) 



n 



n 

7^(P2 -Pi) 
n 



+ l)+cos 



7r(pi + ga) 



+ cos 



7r(Pl + gl) 



+ p 



-1/2 



(9.87) 
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